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^— ^ ■ Cosmological evolution of the QSO luminosity functions (LFs) xat NIR/optical/X-ray bands for 

' 1.3<z<3.5 is investigated based on the realistic QSO spectra. The accretion-disk theory predicts that 

although QSO luminosities only depend on mass-accretion rate, AI, QSO spectra have a dependence on 
black-hole mass, Mbh, as well. The smaller Mbh is and/or the larger M is, the harder becomes the QSO 
[~~. ' NIR/optical/UV spectrum. We model disk spectra which can reproduce these features and calculated LFs 

for redshift z ~ 3 with the assumption of new-born QSOs being shining at the Eddington luminosity. The 
main results are: (i) the observed LFs at optical and X-rays can be simultaneously reproduced, (ii) LFs 
at optical and X-ray bands are not sensitive to Mbh 7 while LFs at NIR bands are; about one order of 
ly-v ' magnitude difference is expected in volume number densities at i/.j ^ lO^^ergs""'^ between the case that 

ly^ I all QSOs would have the same spectral shape as that of A/bh = IO'^Mq and the case with Mbh = IO^^Mq. 

CN ■ (iii) The resultant LFs at NIR are dominated by IO'^Mq black-holes at L/,j < 10^''ergs"\ and by lO^M© 

C^ I black-holes at Lij > lO'^^ergs^^. Future infrared observations from space(e.g.NGST) will probe cosmolog- 

ical evolution of black hole masses. For redshift 2: < 3, on the other hand, the observed optical/X-ray LFs 
can be fitted, if the initial QSO luminosity Lq is below the Eddington luminosity L^dd- Interestingly, the 
best fitting values of -^ = io/-^Edd are different in B- and X-ray bands; £b ~ 2.5 1^- The reason for this 
discrepancy is briefly discussed. 

Key words: accretion, accretion disks — black holes — cosmology: theory — galaxies: active — 
2 '. galaxies: Seyfert 

1. Introduction to note similar evolutionary behavior of the cosmic star- 

- ^^ . formation history (Madau et al. 1996; Connolly et al. 

X ■ 1997; Glazebrook et al. 1999; Steidel et al. 1999), indi- 

H : The cosmological evolution of QSOs and QSO black mating a close link between the QSO evolution and the 

. . . , holes attracts many researchers recently thanks to the cosmic star formation history (Franceschini et al. 1999). 

rapid progress in observational studies of distant QSOs, 

which promotes intensive discussion on their forma- LFs of bright QSOs are well approximated by double 

tion mechanism and formation epochs (e.g., Rees 1984; power laws (e.g. Boyle et al. 1988). This shape con- 

Haiman & Loeb 1999: hearafter IIL99; Kauffmann & trasts those of normal galaxies which show exponential 

Haehnelt 2000; Nulson & Fabian 2000). The efforts have decline at high luminosities. Pioneered by Cavaliere & 

been concentrated on understanding the specific evolu- Szalay (1986) and Efstathiou & Rees (1988), theoretical 

tionary behavior of the bright QSO population. It shows modeling of the QSO LFs has been attempted by several 

a rapid increase from the present day back to redshift authors; e.g., Haehnelt et al. (1998, hereafter HNR98), 

z ~ 2.5 (Schmidt, Green 1983). It is still controversial Haiman & Loeb (1998, hereafter IIL98; IIL99), Cava- 

whether the comoving density of QSOs remained high at here & Vittorini (1998, 2000), and Kauffman & Haehnelt 

z > 2.5, as is suggested by ROSAT X-ray studies (Miyaji (2000). We specifically pick up models by HNR98 and 

et al. 1998, 2000) or rapidly decreases toward higher z, HL98, who consider formation of black holes in hierar- 

as is obtained by optical and radio survey(Pei 1995, see chically growing dark halos, whose formation rate is as- 

also Osmer 1982; Shaver et al. 1996,1999). This apparent sumed to be described by the Press-Schechter theory. It 

discrepancy may simply imply that X-ray surveys could is also assumed that each dark halo finally produces one 

detect optically faint QSOs whose population shows no supermassive black hole with mass being a function of 

rapid decline beyond z ~ 2.5 (HL99). It is interesting mass of its host dark halo, and that QSO black holes 
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shine at Eddington luminosity for certain period by ac- 
creting environmental gas and then quickly fade out as 
accretion material is depleted. The basic shapes of the 
LFs in either B-bands or X-rays are then nicely repro- 
duced. 

However, these studies assume rather ad hoc accretion 
spectra which require improvement by using realistic ac- 
cretion flow models. This is the motivation of the present 
study. We, here, especially focus our discussion on how to 
extract information regarding the cosmological evolution 
of black-hole mass from multi-wavelength LFs. 

The basic accretion theory tells that black-hole accre- 
tion produces energy output of the order of O.lMc^ (with 
c being speed of light); that is, the disk luminosity, L, 
does not contain information regarding black-hole mass 
except for the lower limit of A/bh set by the constraint 
that L should not largely exceed the Eddington lumi- 
nosity, LEdd = 1-5 X 1O3®(MBH/M0)erg s"^ We wish 
to emphasize, however, that spectral energy distribution 
certainly depends on black- hole mass. This is observa- 
tionally clear, since galactic black-hole candidates with 
A/bh ^ IOM0 exhibit distinct spectra in optical to soft 
X-ray ranges from those of typical AGNs. Further, there 
is growing evidence that narrow-line Seyfert 1 galaxies 
(NLSls), which are believed to harbor relatively less mas- 
sive black holes with A/bh — 10^ — IO^A/q, show unique 
soft X-ray features; i.e., enhanced soft X-ray excess and 
large spectral index in X-rays (Pounds et al. 1996; Otani 
et al. 1996; Boiler et al. 1996; Brandt et al. 1997). 

These features can be basically accounted for in terms 
of the standard-disk theory (Shakura & Sunyaev 1973) 
and that with some extension (Mineshige et al. 2000). It 
is well known that the effective (or surface) temperature 
of optically thick (standard) disks depends on ATbh as 
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that is, the peak frequency of the disk spectral peak shifts 
in proportion to Afgjj' x (L/^Edd)^^"'- The smaller A/bh 
is and/or the larger M (or L) is, the harder becomes the 
QSO optical/UV spectrum. X-ray spectra, in contrast, 
practically have no mass dependence, since the electron 
temperature of optically thin, hot plasmas near black 
holes is in any cases kept around 5 x 10^ K regardless 
of black hole mass. This is explicitly demonstrated in 
the framework of a model of optically thin, advection- 
dominated accretion flow (ADAF; Ichimaru 1977, see 
Manmoto et al. 1997 for the spectrum without self- 
similar assumptions). 

Then, we may be able to extract information concern- 
ing black- hole masses as a function of redshift, z, through 
the comparisons of QSO LFs at two (or more) differ- 
ent wavebands; e.g. X-rays, optical, and infrared bands. 
This will be a final goal of the present study. The plan of 



the present paper is as follows: in section 2 we explain our 
models, which are basically the same as those of HL98 
and HNR98 except for adopting realistic QSO spectral 
models. We then give our resultant LFs at B-, X-ray, 
and NIR bands for z ^ 3 and 2; < 2 in sections 3 and 4, 
respectively. The final section is devoted to discussion. 

2. Our Models 

2.1. Basic assumptions 

It is generally believed that most of the mass in the 
universe is in the form of "dark matter (DM)." Accord- 
ingly, the spatial distribution of luminous objects, such 
as galaxies, QSOs, and so on, is likely to follow that of 
DM. Therefore, it is reasonable to use the volume num- 
ber density of dark halos to investigate that of luminous 
objects, including QSOs. We here adopt two simple mod- 
els for calculating LFs, following IINR98 and HL98. The 
basic assumptions made in both models are as follows: 

2.1.1. Relation between halo mass and black-hole mass 
We assume that each dark halo necessarily has only 
one supermassive black hole at the center, though it is 
still controversial obscrvationally. Then, we have 



A//BH = F(Afhalo) 



(2) 



where Afhaio is the host dark-halo mass and the function 
form F{M) depends on models, as prescribed below: 

In model A (HNR98), the relation between A/bh and 
A/haio is nonlinear; 



Af, 



BH 



Cvl^,, = C'Ml'l{l + zfl^ 
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where C" and C are constants and C is treated as a free 
parameter. The physical meaning of this relation is de- 
scribed in Silk & Rees (1998). Equation (g) gives a criti- 
cal black hole mass to bound gas in its host halo against 
the outward wind from QSO. If black hole mass exceeds 
this upper limit, QSO could expel all the gas from its 
host galaxy, and the black hole will never grow further. 
In Model B (HL98), in contrast, a linear relation be- 
tween A/bh and A/haio is assumed; 



AfBH = eAfh 



alo- 



(4) 



Here, e is a free parameter. This expression is based on 
the observed linear relation among the bulge luminosity 
and the black-hole mass (Magorrian et al. 1998, hereafter 
MG98; see also Laor 1998; Ferrarese & Merritt 2000; 
Gebhardt et al. 2000). 

2.1.2. QSO light curves 

We, next, assume that time evolution of the QSO lu- 
minosity at high z ( > 3) follows 



Lit) = LEddGxp 
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where we set t — when a halo collapses. Such a simple 
prescription for a single QSO light curve is known to re- 
produce well the observed LFs at z > 3 (HL98, HNR98). 
If one directly uses the models described above to cal- 
culate QSO LF at z < 2, we over-predict number density 
of QSOs compared to the observed co-moving density of 
QSOs that exhibits rapid decline towards z — from 
z ~ 2.5. Therefore, model assumptions should be re- 
examined. The rapid decay in the QSO density seems 
to be caused by the depletion of fueling mass; i.e., envi- 
ronmental gas surrounding a black hole is not enough to 
shine the black hole at LEdd- Hence, it is reasonable to 
assume that the QSO luminosity should be less than L^dd 
even at the epoch when its host halo collapses (Mcleod 
et al. 1999; Haiman & Menou 2000). We thus introduce 
a parameter, £ = L^/LEdd < 1, and, assuming a relation, 
Lq oc i^Edd oc A/bh, we set 



L{t) = Loexp - 



t 



£Mbh 9{t). 



(6) 



Considering the evolution of gas content in galaxies due 
to star formation activity, (. will be decreasing towards 
z-0. 

2.1.3. Formation rate of dark halo 

In previous studies, the formation rate of dark halos 
and black holes was often regarded as the time derivative 
of Press-Schechter (PS) mass function (Press & Schechter 
1974); 

dAlBudz e dz dMhaio 

At later times (especially at z< 2), however, this is in- 
adequate, since equation (0) gives a negative formation 
rate for small halo masses. This is because that the time 
derivative of PS mass function contain two distinct terms: 
i) the formation rate (which is positive) of halos with 
mass M from objects with lower masses and ii) the de- 
struction rate (which is negative) due to merging of halos 
with mass M into objects with higher masses, and the 
latter dominates over the former at later times for small 
halo masses. Many authors have attempted to derive the 
genuine formation rate in various ways (e.g., Lacey & 
Cole 1993, Sasaki 1994). We here adopt a model pro- 
posed by Kitayama & Suto (1996). Then, we regard the 
formation rate of black holes to be 

d^NBH 1 d dNi, 



dMBudz e dz dM]^, 
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xp{z',z). 
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where, d{dN{orm/dMiiido)/dz is the formation rate of dark 
halos, and p{z' , z) is the survival probability, probability 
that the dark halo formed at z' remains at z without 
merging into objects of higher masses. Kitayama & Suto 
(1996) calculated the genuine formation rate by the merg- 
ing rate of halos of < Af/2 into a halo of mass M. 



After all, the number of free parameters is three for 
each model: QSO lifetime, tq, the initial Eddington ra- 
tio, i, and the constant related to Mbh, C (model A) 
or e (model B). These are to be determined such that 
the model should reproduce the observational data most 
successfully. 

2.2. QSO Luminosity Functions 

Here, we let $(L, z)dL be the number of QSOs per unit 
comoving volume at redshift z, whose absolute bolomet- 
ric luminosity is between L and L + dL. On the basis 
of these assumptions described in subsection 2.1.1-2.1.3, 
QSO LFs, ^{L,z)dL, can be calculated by the summa- 
tion of luminosities of all the QSOs whose luminosity is 
L at redshift z; that is 



$(L,z) 



dMvividz' 
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where 6 is delta-function, g(t) is defined in equation (0), 
and tz^z' is the time between the epochs of redshifts z and 
z'. Here, the lower hmit of Afhaio is Afhaio > 10^Mq[ {1 + 
z)/10 ]^^'^. This corresponds to the virial temperature, 
Tvir > WK. This is the condition for the gas to sink to 
the black hole by cooling effect. 

Practically, it is convenient to define £band, the obser- 
vational band luminosity. We should note that iband is 
a function of Mbh and M in accordance with mass and 
mass-flow rate dependences of the QSO spectra. We, 
hereafter, use '^'(ibandi^) rather than <I>(i,z) and sim- 
ply write ^{Lhand, z) as $(L,z) below. 

The LFs certainly depend on the adopted cosmolog- 
ical model. Though it still remains an open question 
whether the universe is open or closed, recent observa- 
tions seem to indicate a flat universe (e.g., de Bernardis 
et al. 2000). Thus, here we adopt a ACDM cosmology 
with tilted power spectrum, {no,flA,ni^,h,agi^-i,n) = 
(0.35, 0.65, 0.04, 0.65, 0.87, 0.96), as is suggested by of Os- 
triker & Steinhardt (1995). 

2.3. Adopted QSO Spectra 

As discussed before, we need to give disk spectra to 
calculate band luminosities, iband, and thus QSO LFs. 
Here, we adopt a simple but reasonably realistic model 
of QSO spectra as displayed in figure [l] These spectra 
are the sum of three components; emission from a disk- 
corona structure (at A < 1/ini or log;^ > 14.5), that from 
an outer optically thick disk (around A ~ 1/ini), and the 
IR bump (at A > 1/xm). In addition, we simply assume 
strong photo-electric absorption due to IGM (Especially, 
Lyman limit systems mainly contributes to absorption 
here. (Stengler-Larrea et al. 1995)) in the range, 60 
A< A < 912 A. 
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Fig. 1.. The adopted QSO spectra. The sohd lines are the 
cases for M = 0.5 X X/Bdd/c^i while the dotted lines are for 
M = 5.0 X I/Bdd/c'^. Three lines of each set represent the spec- 
tra for Mbh = 10'^\ 10^, IO'^Mq from the top downwards. 



2.3.1. Disk- corona component 

The disk-corona model spectrum is the main compo- 
nent at A < Ifim. Here, we use a model by Kawaguchi et 
al. (2000), since they could, for the first time, reproduce 
the observed broad-band (optical to hard X-ray) spec- 
tral energy distributions of QSOs. The idea is to cou- 
ple a standard-type disk body with Tcff ^ 10^(r/rg)^'^/'' 
K and an advection-dominated disk corona, in which 
Teicc ^ 10^ K. Following Haardt & Maraschi (1991), they 
assume that a fraction, 1 — /, of energy is dissipated in 
the disk body, while the remaining fraction, /, is in the 
coronae. By solving the hydrostatic balance and ID ra- 
diative transfer including inverse-Compton processes for 
each radius, they finally obtain the spectrum at 4 — SOOrg. 
According to this model, the big blue bump (BBB) is by 
thermal black body emission from the disk body at small 
radii, the soft X-ray excess is inverse-Compton scattering 
of the BBB soft photons in the corona at small radii, and 
the hard X-rays are bremsstrahlung radiation from the 
coronae at large radii. 

Free parameters involved with this model are Mbh: 
M, /, ttci Tc, tq, Rin, Rout] black-hole mass, accretion 
rate, fraction of energy dissipated in the corona, viscos- 
ity parameter in the corona, coronal optical depth, disk 
optical depth, inner and outer edges of the corona, re- 
spectively. The calibration of this model spectra is made 
so as to reproduce the so-called QSO composite spectrum 
(Zheng et al. 1997; Laor et al. 1997), in which the spec- 



tral indices (L^ (x v ") are a ^ 0.3 at optical (1 fim 
> A > 2500 A), a ~ 1.0 at UV(2500 A> A > 1000 A), 
a - 1.8 at FUV(1000 A> A), a - 1.6 at soft X-ray(0.2- 
2.0 keV), and a -- 0.7 at hard X-ray (> 2.0 keV), and cut- 
off energy for hard X-ray is about 100 keV. Kawaguchi 
et al. find that one can reasonably reproduce the ob- 
servational QSO composite spectrum with parameters of 
(/, ac, Tc, TO, Rin, Rout) = (0.6, 1.0, 0.6, 1000, 4.0i?g, 
300i?g), Mbh = SxIO^Mq and Af = O.SxLedd/c^ Note 
that M of 12LEdd/c^ corresponds to a disk with £Edd- 
The frequency at which the BBB reaches its peak flux 
varies as Mgjj as in the standard disk model, since the 
main component of BBB is the thermal emission from 
the optically thick disk [see equation (|l|)]. 

Mass accretion rate is assumed to change according 
to equation ^. Thus, the spectral shape also changes 
with the time elapsed after a halo collapses. The spec- 
tra for different A/bh and M are calculated according to 
the disk corona model with the same parameter set given 
above. As M decreases and/or Mbh increases, the spec- 
tral peak shifts towards longer wavelength regimes, as 

A//iEdd and A^bh (Figure 1). These dependences 

are the same as that of the standard disk. Since opti- 
cal/NIR spectral index a is ~ 0.3, luminosity those bands 
varies as Afg^ with fixed M/L-^dd- The X-ray luminosity 
is in proportion to A/bh, on the other hand. 

2.3.2. Outer optically thick part 

In addition, we consider the spectrum of the standard 
disk at larger radii (r > 300 Rg), assuming that disk 
corona only exists inside 300 -Rg. This component ap- 
pears between the IR bump and the BBB in the spec- 
trum. For the total disk size of 1000i?g, the spectral in- 
dex of the standard disk is a '^ —0.3 at this wavelength. 

2.3.3. IR bump 

The second component is the IR bump, typically ex- 
tending from (20-30) to 1 fj.m (Telesco et al. 1984; 
Radovich et al. 1999). At the present, it is widely be- 
lieved that the dust thermal emission is responsible for 
this bump. A supporting evidence is that the higher en- 
ergy end of the bump is found at ~ 1 /im, corresponding 
to a temperature of ~ 1800 K, at which dust sublimates 
(e.g., Kobayashi et al. 1993). Here, we assume that the 
spectral indices of a = 1.4 at A < 10 ^m (Neugebauer et 
al. 1987; PoUetta & Courvoisier 1999), and a = -2.0 at 
A > 10 /im (Rayleigh- Jeans regime), and that the total 
IR power varies with Af , keeping the luminosity ratio of 
the BBB to the IR bump. Although this assumption is 
rather uncertain, it is worth noting that this component 
contributes not to band luminosity, but only to bolomet- 
ric luminosity in the present study. 
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Fig. 2.. The fraction of the luminosity in each band among 
bolomctric luminosity, /bandi ^^ ^ function of Mbh at 2 = 3.0 
and M = 5.0 X Lg^jj/c-^. Line represent B-band(solid hne), 
X-ray(dotted line), I-band(short dash line), and J-band(long 
dash Une), respectively. 



3. QSO Luminosity Functions at z ^ 3.0 

Once the spectra are given, we can calculate the frac- 
tion of the luminosity at a certain band among the bolo- 
mctric one, /band = L\ja,nd/L. The previous models adopt 
the same spectral shape taken from the averaged observa- 
tional one (e.g., Elvis et al. 1994) for wide ranges of Mbh 
and M ; that is, the value of /band was kept constant. This 
treatment is, however, inappropriate, since QSO spectra 
should have some A/bh and M dependences as we saw 
in the previous subsection. In other words, /band should 
change as a black hole grows and M decays with the time. 

Figure ^ summarizes the resultant /band in B-(4035 - 
4765A), l-(0.94 - 1.12 /xm), J-bands(1.15 - 1.33 /im), and 
X-ray (0.5 - 2.0 keV, corresponding to the ROSAT en- 
ergy band). This figure shows that in the 1- and J-bands 
/band dccreascs by one order of magnitude as Mbh de- 
creases by four orders, whereas in the B-band and X-ray 
/band varics less. This distinct behavior in the infrared 
bands compared with other bands arises from the fact 
that EBB spectral shape moves towards shorter wave- 
length with decreasing Mbh Accordingly, this produces 
interesting features in the shape of the LFs. 

The theoretical QSO LFs can be calculated accord- 
ing to equation (^ for each waveband. As noted be- 
fore, the parameters were determined such that the the- 
oretical models should reproduce the observational, rest- 
waveband luminosity functions in B-band (Pei 1995) and 



Fig. 3.. QSO Luminosity Functions calculated based on model 
A. Here, model parameters are iq = 4.0 X lO'^ yr, £ = 1, and 
C = 4.0 X 10^. Note that z = 3.5 for X-ray ROSAT band, 
and otherwise z = 3.0. Top left; X-ray ROSAT band. The 
filled squares are data by Miyaji et al. (1998) for the redshifts 
of 2.3 < z < 4.6. Top right; B-band. The filled squares are 
data by Pei (1995) for the redshifts of 2.5 < z < 3.5. Bottom 
left; I-band. Bottom right; J-band. In each figure, dotted lines 
represent the case that QSOs have the same spectrum (and thus 
the same /band) cis that for black hole mass of A^bh — 10^ and 
lO^^M© at M = 5.0 X I/Edd/c^ (see Figure 1), and the solid 
lines represent the case that each QSO has MBH-dependent 
spectrum. 



in X-rays (Miyaji et al 1998). In the present study, we 
assume that the AGN obscuration play only a minor role 
in the high-z (> 1), high luminosity {L^ > lO"''*'^ erg 
s^^) regimes for the following reasons. It is indeed true 
that there are several intrinsically obscured QSOs ob- 
served (obscured by some material, say, dusty tori and/or 
starburst regions surrounding the active nuclei) . Brandt 
(1997) noted the existence of type 2 QSOs and Veilleux 
et al. (1999) pointed out that a part of ULIRG is the can- 
didate for highly obscured QSOs. However, it is unlikely 
that they contribute significantly to the QSO LFs at the 
bright end, since Akiyama et al. (2000) reported that 
type 2 AGNs known to date are dominated by nearby 
low-luminous objects in 2 - 7 keV band (see their Fig. 9). 
Further, they plot the redshift distributions of type 1/2 
AGNs as functions of redshifts (see their Fig. 10b). Then, 
they have calculated the expected redshift distribution 
of type 1 AGNs using the hard-band X-ray LFs of type 
1 AGNs (Boyle et al.l998), finding that this expected 
distribution agrees with that of their samples of type 1 
AGNs (obtained by ASCA Large Sky Survey) by the pos- 
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Fig. 4. 



Same as figure pJ but for model B. Model parameters are 
= 3.5 X 10^ yr, £ = 1, and e = 4.0 X 10""'. 



sibility of 64%. Next, they did the same but for type 2 
AGNs assuming that the shape of LFs of type 2 AGNs 
is the same as that of type 1 AGNs. The expected num- 
ber of the type 2 AGNs at z > 0.4 is about 10 among 
30 identified AGNs in total, whereas there was no type 
2 AGNs detection in their data at high redshift. That 
is, it is very unlikely that type 2 AGNs are distributed 
in the same way as that of type 1 AGNs; such a proba- 
bility is only 5%. Certainly, there exists a deficiency of 
type 2 AGNs at z > 0.4, at least, in the high luminosity 
regimes. Still, we admit that the source number is in- 
sufficient to derive solid conclusion from their data. We 
should await further observational study to be conducted 
in near future. 

In the present study, we thus consider only type 1 
QSOs in the calculations of the QSO LFs to be com- 
pared with the observed QSO LFs. It is still possible 
that, depending on X-ray energy ranges, the fraction of 
type 2 AGNs may vary to some degree, giving rise to 
slightly different behavior in LFs at soft and hard X-ray 
bands. This issue will be discussed in a future paper. 

The observed luminosity functions depend on the cos- 
mological parameters (Hubble parameter, h, and the de- 
celeration parameter, go), as well as the spectral index at 
the relevant bands (for iiT-correction) . The dependence 
on h and go is simple; '^{L, z) ex dV~^d~j^ , L oc d^ (here, 
dV is the comoving volume element, and d^ is the lumi- 
nosity distance). As shown by Hartwik & Schade (1990), 
the dependence on a is comparatively large. Pei (1995) 
presents data only for {h,qo,a) = (0.5,0.5,0.5) and 



(h,qo,a) — (0.5,0.1,1.0). In our spectral model, a is 
nearly 1.0 in the B-band at redshift, z ^ 3. Thus, we 
modify Pel's data for {h,qo,a) — (0.5,0.1, 1.0) to satisfy 
ACDM. Since the observed optical and X-ray LFs are 
corrected with an assumed a of unity, which is equiva- 
lent to no K-coUection, we calculate theoretical LFs at 
observer's frame. Once the parameters are determined 
by optical and X-ray LFs, then the luminosity function 
expected in the I- and J-band can be obtained. 

Figure 3 and 4 represent the calculated luminosity 
functions for model A and model B, respectively. In all 
the panels of both figures, the dotted lines represent mod- 
els assuming that all QSOs have the same spectrum as 
that with Mbh of 10^ (lower line) or IO^^Mq (upper line) 
and M of 5.0 x L^dd/c^, whereas the solid lines represent 
models that each QSO exhibits different spectra accord- 
ing to variations of black- hole mass. 

As figures 3 and 4 show, the dotted lines nearly coin- 
cide with each other in X-ray and B-band, whereas they 
can be clearly separated in I- and J-band. In other words, 
LFs in X- and B-bands contain practically no informa- 
tion regarding black-hole mass, while LFs in NIR bands 
do. This can be explained in terms of the different be- 
havior of /band ill those band, as are illustrated in figure 
2. In the X- and B-band, /band varies little (within a 
factor of ^2), even if A/bh changes by four orders of 
magnitude, while in I- and J- bands /band appreciably 
decreases (a factor of ^ 10) as Mbh decrease. In the 
lower two panels of figures 3 and 4, especially, the solid 
line approaches the dotted line of Mbh = IO^^Mq on the 
higher luminosity side (Lband ^ 10^^ erg s~^), whereas it 
approaches that of Mbh — IO^Mq on the lower luminos- 
ity side (iband^lO^'* erg s~^). This clearly points that 
the more massive the black hole is, the brighter the QSO 
becomes at NIR bands (Figure 2). 

The dotted lines in figures 3 and 4 are inconsistent with 
our model spectra, since those lines are plotted under 
the assumption that QSOs with different Mbh have the 
same spectral shape, while we expect MnH-dependence 
of the QSO spectrum as is demonstrated in figure ^ As 
mentioned before, the previous studies have used the one 
observational spectrum for all QSOs. In this sense, the 
dotted lines also represent the cases, in which the spec- 
tral shape is fixed as in the previous models (e.g. IIL98, 
HNR98). The observational mean spectrum by Elvis 
et al (1994) is nearly equal to our spectrum model of 
Mbh - IO^Mq and M = 5.0 x L^dd/c^- The upper 
panels of figures 3 and 4 show that whichever model one 
may use, the previous model or our model, the luminosity 
functions hardly differ in B-band and X-ray. The lower 
panels, on the other hand, demonstrate that black- hole 
mass dependence of the QSO spectra manifests in the 
luminosity function in I- and J-band. Therefore, if QSO 
luminosity functions are observed at z ~ 3 in these wave 
bands in future, we can, in principle, constrain the Mbh- 
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dependence of QSO spectrum, thus probing cosmological 
evolution of the black-hole mass in QSOs. 

Despite the fact that QSO diminishing timescale is 
longer in model A {Iq ~ 4.0 x 10^ yr) than in model 
B (3.5 X 10^ yr), these figures look very similar. This is 
because the mass fractions are different: e — 4.0 x 10"'' in 
model B, while in model A e is not constant but around 
~ 10~^. Here, we note that if we use larger tg and 
smaller e in model B, we cannot reproduce observed LFs. 
Since the inclination of LFs mainly depends on e, theo- 
retical LFs become more steeper with smaller e. 

Which value is more reasonable for e? The ratio of 
the black-hole mass to the host halo mass (e) is not 
directly observable, but MG98 suggest for the nearby 
galaxies and QSOs, the ratio of the central black-hole 
mass to the bulge mass is MEnZ-Mbuigo ^ 0.006. Then, 
the value of e can be estimated through the baryon mass 
fraction in the total, fib/^o — 0.1. If most of baryonic 
matter is contained in the bulge, we have e ~ 10^^'°, 
more consistent with the Model B assumption. MG98 
also support the linear Mbh — -^^haio relation, equation 
(^. On the other hand, Fukugita et al. (1998) noted 
that the bulge density is less than the baryon density by 
about one order of magnitude, Slbuigo — O.lilb- If so, we 
have e ^ 6 X lO^'"', in good agreement with the Model 
A assumption. More recently, Meritt & Ferrarese (2001) 
noticed the tight Mbh — c relation (with a being halo 
velocity dispersion); Mbh oc cr'*-^^, which is close to equa- 
tion (I). They also show A//BH/A^buigc ^ 0.0012, which 
is by a factor of ~ 5 smaller than the value obtained by 
MG98 and could rule out model B. The similar conclu- 
sion was derived based on the mass estimate based on 
the reverbration mapping (Laor 1998). 

Another parameter iq should also constrain models. 
By the usual picture of Eddington-limited accretion, 
black hole only grows as Mbh oc exp(i/iEdd)j where 
^Edd ~ 3 X 10^ yr is the e-folding time. As HNR98 men- 
tioned, this supports Model A. 

4. QSO Luminosity Functions at z ~ 2 

Finally, we present the calculated LFs at intermediate 
redshifts. Here, we only display the cases with Model 
A, since Model B produces similar results except for the 
values of tq (see figures 3 and 4). Model parameters are 
iq — 4.0 X 10 ^r and C = 4.0 x 10^, which are the same as 
those in figure 3 (2; '^ 3.0). Then, by setting ^ to be a fit- 
ting parameter, we can well reproduce the observational 
QSO LF at z = 2.0 and z = 1.3, as are demonstrated in 
figure 5. Here, we let l^ and £b be the initial Eddington 
ratios {£) in X-ray and B-band, respectively. The best 
fit models give (4,^b) == (0.24,0.66) at z = 2.0, and 
(4,£b) = (0.10,0.25) at z = 1.3. Infrared LFs are cal- 
culated by using the same £ as those of B-bands and are 
displayed in figure 6. 
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Fig. 5.. The same as figure 3 (model A) but at different red- 
shifts: z = 2.0 (solid lines) and z = 1.3 (dashed lines), respec- 
tively, and in B-bands (top right). Model parameters are the 
same as those in figure 4 except for the initial Eddington ra- 
tios: e^ = 0.24 {z = 2.0) and 0.10(2 = 1.3), and Ib = 0.66 
(left) and 0.25 (right). Also plotted are the observed LFs. The 
filled squares in the top left panel represent the observed LFs in 
X-ray ROSAT bands taken by Miyaji et al. (1998) at redshift 
of 1.6 < z < 2.3 and at 0.8 < 2 < 1.6, while the filled and open 
squares in the top right panel are data from Pei (1995) for the 
redshifts of 1.9 < z < 2.2 and 1.2 < z < 1.4. 



Similar to figures 3 and 4, mass dependence of LFs 
at NIR bands is evident. The more luminous quasars 
are, the more massive become central black holes. As 
mentioned before, we have used the formation rate of 
dark halo by Kitayama & Suto (1996), instead the time 
derivatives of PS mass function. The distinction between 
the two is considerable in these redshifts: For example, 
in X-ray, PS mass function under-estimates the number 
densities of QSOs at Lx < lO'^^ergs"^ by a factor of 10 
when we choose £ to fit LFs at higher Lx side. This is 
because on the lower L^ side, negative destruction rate 
domiates the time derivative of PS mass function. 

Figure 7 plots the time evolution of the initial Edding- 
ton ratios. Our model requires £b being about twice as 
much as €x at both of z = 2.0 and z = 1.3, which con- 
trasts the cases at z ^^ 3, in which the same Eddington 
ratios, .^b = -^x = 1, can give reasonable fit to the obser- 
vations. This implies that for lower mass- flow rates at 
relatively low redshifts X-ray radiation tends to be sup- 
pressed in the observations with respect to our modeled 
value for the same optical radiation. In other words, our 
model spectra may over-predict X-ray fluxes at £ ^ 1, if 
the discrepancy in £b and £^ is real. However, it is too 
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Fig. 6.. QSO LF at I- and J-band at intermediate redshift. Dot- 
ted lines represents the case that QSOs have the same spec- 
trum as that for black hole mass of Mbh — lO^andlO-'^-'^MQ at 
M = 5.0 X i^Edd/c^- Here, we assume the initial Eddington 
ratio to be equal to that at B-band, £ = £b- 



early to take this difference seriously, since the observa- 
tional selection criteria of QSOs are somewhat different 
between X-ray and B-band LFs. Especially, error bars 
in X-ray are comparatively large, since smaller number 
of QSOs were used at X-ray LFs. We should await fu- 
ture intensive observations to be performed (see the next 
section for discussion). 

Several issues still remain to be solved. First, the as- 
sumption of ^ = 1.0 at z ^ 3.0 may not be obvious. 
This is a naive expectation, and high redshift QSO sur- 
vey may point that £ < 1.0 even at z ^ 3. Next, we 
have not considered the time evolution of iq and e. We 
also did not consider depletion of fuel mass (Kauffman, 
Haehnelt 2000; Cavaliere, Vittorini 2000). These evolu- 
tion may explain decreasing LFs at z < 2 without the 
evolution of £. 



5. Discussion 

We construct a simple model for LFs based on the pre- 
vious models by HL98 and HNR98 but newly adopting 
the realistic disk spectra which have Mbh and M de- 
pendence. Our models for LFs can give reasonable fit to 
the observed LFs both in B-band and X-rays for redshift 
z ~ 3. We also show that these LFs are not sensitive to 
black-hole mass distribution in QSOs, but instead LFs 
at NIR, which will be available in the near future by IR 



Fig. 7.. The time evolution of the best fitting initial Eddington 
ratios, <?b (dotted line) and <?x (solid line). 



QSO surveys, should contain information as to Mbh- 

Our model could also reproduce the shape of LFs at a 
variety of redshifts less than z = 3.0, however, we need 
to assume sub-Eddington initial luminosity of black holes 
at their formation epoch. This is consistent as the results 
obtained by Haiman & Menou (2000) who considered the 
B-band LFs only. From the semi-analytical approach for 
galaxy/QSO formation, Kauffmann & Haehnelt (2000) 
suggested that the sharp decline in QSO number densi- 
ties from z ~ 2 to z ~ results from a combination of 
three factors: (i) a decrease in the galaxy-galaxy merging 
rate, (ii) a consumption of gas due to star formation (cor- 
responding to decreasing £), and (iii) an adopted assump- 
tion of increasing accretion timescale (i.e., decreasing £ 
and increasing tq). 

Our model predicts ^x < ^b by a factor of ~ 2.5. This 
may indicate that our model spectra overestimate X-ray 
flux at luminosity L <C iEdd- In other words, we need 
to suppress X-ray flux compared with B-band one in the 
model spectra at low luminosities. This may require sig- 
nificant improvements of the model spectra adopted in 
the present study. Cosmological evolution of QSO spec- 
trum requires further research. We also note that this 
may take part in the discrepancy of observed comoving 
density of QSOs at z > 2.5 between B-band and X-rays. 
If QSO accretion-disks turn into the ADAF regime at 
low redshift due to low accretion rate (e.g., Yi 1996; 
Haiman & Menou 2000), optical radiation should be 
much suppressed compared with X-ray radiation; the 
broad-band spectrum of ADAF is generally harder than 
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the present disk-corona spectra, on the contrary to our re- 
sults. Intrinsic dust extinction at (rest) optical/UV leads 
to the reduction in observed B-band luminosity, which 
again produces the opposite effect. Rather, our result 
may suggest a drastic shrinkage of the X-ray emitting re- 
gion (probably accretion disk corona) at low luminosities. 
To probe this, we need more studies concerning the con- 
ditions creating and/or sustaining hot coronal flow above 
cool accretion disk. 

In the present study, we calculated the formation rates 
of dark halos by the formalism based on PS mass function 
(PS halo), assuming each PS halo harbors only one black 
hole. Generally speaking, however, PS halos are not pre- 
cisely equal to the galactic halos. For example, the mass 
scale of Mhaio ^ 10^'^Mq corresponds to a galaxy cluster 
and, hence, it is not correct to assume that such a huge 
PS halo contains only one black hole (this is known as 
the over-merging problem). Thus, we should be careful 
with the maximum halo mass Mhaio used in the calcu- 
lation of QSO LFs so as not to exceed a certain limit. 
At z — 3.0, the halos which mostly contribute to LFs 
at the bright end (for B-band, Lb ~ lO^'^'^iB.sun) have 
A^haio ~ IO^^Mq, but at z = 1.3 the maximum halo 
mass attains 10^^ Mq. Therefore, our calculated LFs suf- 
fer over-merging problems only at higher luminosity side 
at lower redshift. But to evaluate £, we can conclude that 
this effect is small. 

Throughout the present study, we assume that each 
dark halo contains only one massive black hole whose 
mass is scaled with the host dark-halo mass. It is inter- 
esting to note recent discovery of multiple, intermediate- 
mass black holes being created in post-starburst regions 
in M82 (Matsumoto & Tsuru 1999; Matsushita et al. 
2000). This is not consistent with our assumption, but 
it may be that such intermediate-mass black holes will 
eventually merge into one big hole within timescales 
much shorter than the Hubble time. 
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